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Abstract 

Assuming that there exists a species of heavy sterile neutrinos (m u > 1 
keV) and their decays can serve as a heating source of the hot gas in galaxy 
clusters, we study how the observational constraints on cooling flow limit the 
mass of these sterile neutrinos. We predict a relation among the luminosity, 
total mass and the redshift of a cluster, and we compare this relation with 
data from 12 clusters to obtain an estimate of the decay rate of the sterile 
neutrino. 
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1 



Typeset using REVT^X 



I. INTRODUCTION 



Recent observations indicate that there are not enough cooling flows in galaxy clusters 
to explain why hot gas in some clusters are so hot [1], [2]. On the other hand, neutrino 
oscillation experiments (solar, atmospheric and LSND experiments) indicate that 3 ordinary 
neutrinos cannot explain the three different scales of the mass squared differences between 
neutrino species Am 2 's, of about 10~ 5 eV 2 , 10~ 2 eV 2 and 1 eV 2 [3], suggesting that one or 
more sterile neutrinos must be added. In this article, we study the possibility that a species 
of massive sterile neutrinos (m„ > 1 keV) exists and their decays provide a heating source to 
maintain the temperature of the hot gas. We use this model to predict a relation among the 
luminosity, total mass and redshift of a cluster, and we obtain an estimate of the decayrate 
of the sterile neutrinos by using available data from 12 clusters. This idea of the radiative 
decay of neutrinos as an energy source has been proposed to account for reionization [4]. 

In standard cosmology, neutrinos with m u < 1 MeV are classified as hot dark matter 
(HDM) because they are relativistic at decoupling. As the universe expands, they cool down 
and become non-relativistic after a redshift of 



where T u0 is the neutrino temperature today. In the non-relativistic regime, the average 
neutrino speed at neutrino temperature T u is [5] 



They do not contribute to structure formation if their average speed is too high. For example, 
at matter-radiation equality (redshift ~ 10 4 /i 2 , h being the Hubble parameter), the neutrino 
speed is close to c if m v <§C 10 eV. By using perturbation theory on neutrino clustering, it 
was found that massive neutrinos can gravitationally cluster if m v is greater than several 
eV's [5]. Therefore, neutrinos with mass m v >10 eV can be classified as warm dark matter 
(WDM) rather than HDM. 




(1) 




(2) 
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II. STERILE NEUTRINO DECAY IN CLUSTERS 



In this model, we assume that all 3 types of active neutrinos with masses m 1? m 2 and m 3 
did not contribute to structure formation and are thus classified as HDM [6]: 

mi + m 2 + m 3 2 

94 eV = Q " h ' (3) 

where VL V is the light neutrino density parameter. In addition, we assume that all cold dark 
matter are composed of two types of sterile neutrinos with masses m sl and m s2 , and the 
latter will decay into mi neutrinos with decay rate T: 

v S 2 v\ + 7- (4) 

The massive sterile neutrinos may be produced in the early universe by oscillations from 
active neutrinos [7] . From the standard cosmology, the number densities of sterile neutrinos 
n s \ and n s2 at present (time = to) are upper bounded by [4] 

n , 1 = L2xl0 -* n Jh^)(^)(") = ^L no , (5) 



m al I V0.23/ \ 0.72 / m sl 



where a± < m s \ and a 2 < m S 2 are constants which are dependent on the cosmological 
sterile neutrino density parameters Q s i and Q S 2 respectively, and no is the number density 
of one type of the active neutrinos now (n ~ 100 cm -3 if there were no sterile neutrino 
decays [8]). If the sterile neutrinos can decay to mi neutrinos emitting a photon with energy 
e = (m s2 — mi)c 2 ~ m s2 c 2 , the number density of the mi neutrinos now is given by 



ni = n 



l + _^_(l_ e -rto) 
m s2 



(7) 



Therefore, we can write the cold dark matter density parameter Q m as 



ai + a 2 e rt0 

Qm ~ 94 eV • (8) 
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The interaction of the radiated photons and electrons in the cluster hot gas can be treated 
as a heating source of the hot gas. There exists a lower limit of the energy of the photons 
such that they can undergo Compton Scattering with electrons rather than Inverse Compton 
Scattering (see Appendix). The energy of photons must be greater than the energy of the 
electrons in order to transfer energy to the hot gas, which typically has temperature of about 
several keV's in clusters. Therefore, in our model, m s2 must be at least of keV order. Here, 
we set mi = m 2 = m 3 , as neutrino oscillation experiments indicate that the difference in 
the mass squared of three active neutrino species, Am 2 's, are about 1CT 5 eV and 1CT 2 eV 
[9]. 

Eq. (8) gives one relation between two unknowns, a\ and a 2 . Another relation comes 
from the decay of m s2 neutrinos, which we assume to be the source of hot gas luminosity. 
The total non-baryonic mass in a cluster is given by 

M = Mi + M 2 + M 3 + M sl + M s2 , (9) 

where Mi, M 2 , M 3 , M sl and M s2 are the total masses of three types of active neutrinos and 
two types of sterile neutrinos respectively. By combining Eqs. (5), (6), (7) and (9), we get 

M 



N s2 = — 



a 2 e- rto 



(10) 



m 1 [3r s2 mi + a 2 r s2 e~ rt ° + a 2 (l - e~ r '°) + a\r s2 _ 
where N s2 is the total number of sterile neutrinos inside a cluster, and r s2 = m s2 /mi. The 
total power given out by the radiative neutrinos in a cluster at time t is L x — N s2 Te = 
N s2 Tr s2 mic 2 . Therefore we get 

,-rt 



MIV 



a 2 e 



(11) 



_3m 1 + a 2 e r * + a 2 (l-e rt )/r s2 + «i_ 
At the time of CMB last scattering (ti 10 13 s), the cold dark matter density pi is given 

by 



pia t = ain + a 2 n e \ (12) 

where ai is the scale factor of the universe at ti. In order to heat up the hot gas, the value 
of T must be greater than 3 x 1CT 17 s _1 by using Eq. (11) (also see Fig. 1). Therefore the 



factor e r * ; R3 1 and p^af « ain + a 2 n . At the time of galaxy formation (t g > 10 17 s), the 
dark matter density p g is 

p g a 3 g = «in + a 2 n e~ rts , (13) 

where a 9 is the scale factor of the universe at t g . The factor e~ vtg < 0.03 is much less than 1, 
and we get p g a 3 ~ ain . Therefore, at t g , m s i neutrinos dominate all the dark matter. We 
can therefore obtain another relation between cci and a 2 by taking the difference of Eq. (12) 
and Eq. (13). So the density of m s2 neutrinos p s2 at t\ is given by 

333 

Pfggj = PlQ-i _ Pg a g _ qCMB _ c,2dFGRS ^ a 2^0 /^x 

„ rn m ~ „ > V / 

Pc Pc Pc Pc 

where p c is the critical density of the universe, f2^ MB /i 2 = 0.14 ± 0.02 is from WMAP data 
and n 2 ^ FGRS h 2 = 0.134 ± 0.006 is 2dFGRS data corresponding to CMB last scattering and 
galaxy formation respectively [10]. The corresponding value of a± = 94 x f2^ FGRS /i 2 eV = 
12.60 ± 0.56 eV, which means that m s i must be greater than 12 eV and consistent with our 
assumption that they can form structures. The upper bound of a 2 obtained by Eq. (14) is 
about 3 eV. 

By using these values, we can find the value of Y as well as the corresponding mass 
components in cluster. Since a.\ is much greater than mi, a 2 e~ rt ° and a 2 /r s2 , we can 
rewrite Eq. (11) as 



t 







(15) 



where z is the redshift of a cluster and we have assumed matter dominated expansion and 
r s2 is much greater than 1. By using data of redshifts, luminosities and masses from 12 
clusters (see Table 1), we can estimate V from the slope of hiL x /M vs. (1 + z)~ 3 ^ 2 (see 
Fig. 1). The slope and y-intercept obtained are — 18.23±4.19 and 18.50±3.87, which imply 
T~ (4.6 ±1.1) x 10~ 17 s -1 . 

The corresponding mass components in clusters with i — 1,2, 3, si or s2 are given by 

Mi = ^M, (16) 



where fij is the density parameter of neutrinos with ui{. For a typical cluster with M ~ 
1O 15 M , the present value of M s2 is about 1O 6 M . 

III. TEMPERATURE PROFILE OF HOT GAS IN A CLUSTER 

We can predict the temperature profile of the hot gas in a cluster by using a simple 
energy flow equation. We first look at the mechanism of how the energy from neutrinos 
decays is transferred to the hot gas. 

Since the m s2 neutrinos have mass greater than keV, they must have gravitationally 
collapsed to galactic scale structures [11], [12]. The radius of a hydrostatic 'neutrino star' is 
given by [13] 

p on 7 ( M * 2 V V V m * 2 V 8/3 m 
R = 20 - 7PC [WM^) llkevJ • (17) 

From our estimate, M s2 and m s2 are of order 10 6 solar mass and keV respectively. Therefore, 
the size of a corresponding 'sterile neutrino star' is about 20 pc which may hide deeply inside 
the galactic bulge. 

Suppose a photon travels from r = and collides with electrons in a cluster. The number 
of collisions of a photon and electrons within a radius r is approximately given by the optical 
depth: 

rr dy 1 rr 

x= -Y- = / n e crdr', (18) 



o I Jo 

where I is the mean free path of a photon, and n e and a are the electron number density 
and Compton cross section respectively. In a cluster, the electron number density is given 
by 

Ue = [l + ( r /r c ) 2 ] 3 ^/ 2 ' (19) 

where n c , r c and f3 c are parameters in the cluster beta model [14]. In Fig. 2, we plot the 
number of collisions x vs. r. We can see that x is much less than 1 for all r. Therefore, a 
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photon does not collide frequently with electrons in a cluster. However, inside the Milkyway 
galactic bulge, the number density of electrons is 

Tl! 

Ue = [l + (r/r>W- 8 ' m 
where n' c and r' c are parameters which are equal to 1.6 x 10 8 cm" 3 and 0.34 pc respectively 
[15]. Fig. 3 shows the number of collisions inside the galactic bulge. We can see that x is 
much greater than 1 for small r. Therefore, we can believe that the energy of the photons are 
first absorbed by the electrons inside the galactic bulge and then the energy is transferred 
to the hot gas by conduction inside the cluster. The mean free path for conduction is [16] 

^Hwk)' (w^T • < 21 > 

which is smaller than the length scale of a cluster 1 Mpc) with the temperature of hot 
gas T g ks 10 s K. The energy flux per unit time for conduction is given by 

j(r) = -K c <f, (22) 

where 

K c = lMn e l c k — M . (23) 
\m e J 

We can obtain the temperature profile T g (r) by solving Eq. (22) with j(r) = L(r)/47ir 2 and 
L(r) = J Anr 2 e s {r)dr if we know the energy source distribution e s (r) in a cluster, which is 
approximately proportional to the sterile neutrino distribution. We demonstrate possible 
temperature profiles by using two models. In one we assume that the power sources are 
located at the centre of a cluster with e s {r) = S(r)L/ '4nr 2 , where L is the total luminosity 
of the cluster. In the other one, we assume that the power sources are distributed as 

2/2 

e s = Ae~ r /Tc , where A can be fixed by the total luminosity L. Fig. 4 shows the temperature 
profiles of the two models. We can see that the temperature varies only within 6 percents of 
the central temperature T , which is due to the large value of K c . Therefore, the energy is 
transferred from the sources (decayed m S 2 neutrinos) to the hot gas quite efficiently within 
the Hubble time. 



IV. DISCUSSION AND SUMMARY 



It has been a puzzle why some clusters are still so hot even though their cooling times 
are shorter than the Hubble time. Even though cooling flow models [16] can help to explain 
this phenomenon, significant cooling flow in clusters has not been observed. This is known 
as the 'mass sink problem' [1]. Also, the cooling flow model predicts that the temperature 
of the hot gas should be quite inhomogeneous. However, the observed hot gas temperature 
profiles are quite homogeneous and gas with temperature less than 1 keV does not exist in 
the amount predicted. Donahue and Voit summarized the drawbacks of the cooling flow 
models and suggested that any successful model must explain the mass sink problem and 
the homogeneous temperature profile with positive core temperature gradients extending to 
~ 10 2 kpc [1]. Here, we present an alternative solution of the mass sink problem. We have 
discussed how the radiative decay of sterile neutrinos can heat up the hot gas in clusters. 
Our model predicts a redshift dependence of L x /M, Eq. (15), which is consistent with the 
observed data (see Fig. 1), and we obtain the decay rate of the sterile neutrinos in clusters, 
which is about T ~ (4.6± 1.1) x 1CT 17 s _1 . Our picture does not require any cooling flow, and 
the radiative heating power by neutrino decays results in a quite homogeneous temperature 
profile. Therefore, our model can fit the observed information and solve the cooling flow 
problem for a range of T. 

If we can constrain the values of mi and a\ [6], we can get a better picture and fix Y. 
The radius of a 10 13 solar mass m s i neutrino star is of order 10 kpc which is the same scale 
as a galaxy. Therefore we propose that the galactic dark matter may be composed of m s i 
and m s2 sterile neutrinos. By examining the rotation curve of a galaxy, we can obtain a 
feasible range of m s i. If m\ is below 1 eV, the size of an mi neutrino star is greater than 
1 Mpc, which is the cluster scale. Therefore, the active neutrinos may affect the properties 
of hot gas and clusters [12]. Different neutrinos may correspond to structures in different 
scales. It would be interesting to study whether our model can solve both the dark matter 
problem and the cooling flow problem. 
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V. APPENDIX 



We calculate how the Compton scattering of photons and electrons can heat up the 
hot gas in clusters. In the usual Compton scattering calculation, the electron is at rest. 
However, the temperature of the hot gas in most clusters are at 10 7 — 10 8 K, and so we 
should generalize the situation by transforming to the electron's rest frame [17]. Therefore, 
we have 

e' = 7 e(l-/3cos0), (24) 

e s = 7 4(1 + P cos 9' s ), (25) 

where e and e' are the initial energy of the photon in the lab frame and electron's rest frame 
respectively, and 9 is the angle between the initial photon direction and electron velocity 
in the lab frame. The subscript 's' indicates the scattered energy or angle. The scattered 
power in the rest frame is given by 

P' = caj e'J(p' s )d 3 p' s , (26) 

where a, f and p' are the Compton cross section, distribution function and momentum of 
photons respectively. Since the distribution of the photons is not necessarily isotropic, we 
need to transform back to the lab frame and integrate over all 0' 8 . Making use of d 3 p' = 
d 3 p(l — [3 cos 9), and combining Eqs. (25) and (26), we have 

P' = ccrJ 7 2 (1 - (3cos9 s ) 2 e s f(p s )d 3 Ps . (27) 

After summing up all the momenta of the photon, we have 

P' = ca-f (l + | o °° e s n(e s )de s = ca^ (l + ^ U phs , (28) 

where U phg is the scattered energy density of the photons in the lab frame. Since P' — 
and < dE s >= 7 < dE' s >, we have P' = P, where P is the scattered power in the lab 
frame. In the lab frame, the rate of energy removed from the photon field is 



dE f 
P = —jj- = —ca J en(e)de = —caU p h, 



(29) 



where U p h is the energy density of the photons in the lab frame. Therefore, the net energy 
radiated is the difference of energy scattered and energy removed in lab frame: 



2 f P 2 \ 
P net = co-7 I 1 + — I U phs - caU ph . 



(30) 



If all the photons undergo Thomson scattering with electrons, then e' = e' s in Eq. (30) and 
Uph will be the same as U p h s - However, in general, e' may not be equal to e' s . Suppose we 
have only a particular energy for the photon field, U p h = n ph e. From Eq. (30), we have 



P net = can. 



ph 



Y 1 + 



(e - Ae) - e 



where Ae = e — e s . After simplification, Eq. (31) will become 



P net = cariphrf 



From Eqs. (24) and (25), we have 



> 2 



< Ae' >= e' - e' s = 7AC 



(31) 



(32) 



(33) 



Therefore, the net power radiated in the lab frame is 



Pnet = CaUphl 2 



(34) 



If the net power radiated is a positive value, that means there is some energy removed from 
the electron gas and thus the energy of photon field is increased. This is known as Inverse 
Compton Scattering. However, it is possible that the value of the net power radiated is 
negative, which means that the electrons absorb energy from the photon field. If the photon 
and electron Compton scatter in the rest frame of the electron, then we have 



Ae' 



11 22 
-^-r(l-cos0) = -^(1 -/3cos#) 2 (l -cos0). 



(35) 



We integrate over all the angles to get the averaged energy gain in one scattering: 



10 



< Ae' >= 



m e c z 



Combining Eqs. (32), (33) and (36), we get 



P net = can ph Y 



1 + ^ 



2 



(36) 



(37) 



3 m e c* 

Therefore, the criterion for energy absorption from the photon field is 

Am e c 2 (3 2 

' > 37(l + ?)(l + f)' <38> 

Since < B 2 >= -^j, we can write e in terms of kT. For non-relativistic electron, 3 « 1, 
e > AkT/^j. So, if the photon energy is greater than the kinetic energy of the electron, the 
electron will gain energy and vice versa. 
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TABLES 

TABLE I. X-ray Luminosities (L x ), hot gas masses(M 9 ), total masses (M tota i) and redshifts(z) 
of 12 clusters [1], [2]. 
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Fig.l. ln(L x /M) vs. (1 + z) 3 / 2 for 12 different clusters. The solid line is the best fit 
line with slope -18.23 ± 4.19. 
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Fig. 2. Number of collisions of a photon with electrons vs. distance from the cluster 
centre. We assume r c = 100 kpc, n c = 10~ 3 cm~ 3 and (3 C = 0.7. 
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Fig. 3. Number of collisions of a photon with electrons in Milky Way with r' c = 0.34 pc 
and n' c = 1.6 x 10 8 cm" 3 . 




Fig. 4. The temperature profiles of the two models discussed in the text. The solid line 
represents the model with the sources at the centre of the cluster. The dashed line represents 
the model with e(r) = Ae^l^ . We assume L = 10 44 ergs s -1 , r c = 100 kpc, and the central 
temperature T = 10 8 K. 
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